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Abstract 

If supersymmetry (SUSY) will be discovered, successful models of flavour not only have to 
provide an explanation of the flavour structure of the Standard Model fermions, but also 
of the flavour structure of their scalar superpartners. We discuss aspects of such "SUSY 
flavour" models, towards predicting both flavour structures, in the context of supergravity 
(SUGRA). We point out the importance of carefully taking into account SUSY-specific 
effects, such as 1-loop SUSY threshold corrections and canonical normalization, when 
fitting the model to the data for fermion masses and mixings. This entangles the flavour 
model with the SUSY parameters and leads to interesting predictions for the sparticle 
spectrum. We demonstrate these effects by analyzing an example class of flavour models 
in the framework of an SU (5) Grand Unified Theory with a family symmetry with real 
triplet representations. For flavour violation through the SUSY soft breaking terms, the 
class of models realizes a scheme we refer to as "Trilinear Dominance" , where flavour 
violation effects are dominantly induced by the trilinear terms. 
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1 Introduction 



The flavour puzzle is one of the biggest open questions in particle physics. Currently, one 
may subdivide it into three main parts: the puzzle associated with the masses and mixing 
properties of the quarks and charged leptons, the new puzzles that have been added with the 
discovery of neutrino masses and mixing, as well as the issues connected to CP violation. 

Regarding the masses of quarks and charged leptons it is tantalizing that down-type quarks 
and charged leptons show a similar hierarchical pattern, which differs substantially from the 
pattern of the up-type quark masses. In the lepton sector, in addition to the smallness of 
the neutrino masses the observation of large flavour mixing, in contrast to the small quark 
mixing, has to be explained. Concerning CP violation, there is on the one hand the strong CP 
problem, but on the other hand the observation that the unitarity triangle is, at least, almost 
right-angled, i.e. a ~ 90°. If supersymmetry (SUSY) or any other kind of new physics will be 
discovered at the LHC, this would add another "dimension" to the flavour puzzle. Any kind of 
new physics has to face severe constraints from flavour physics, e.g. flavour changing neutral 
currents in the quark sector (FCNCs) , lepton flavour violation (LFV) as well as electric dipole 
moments (EDMs) . In the case of SUSY this is often referred to as the "SUSY flavour puzzle" . 

Recently a lot of effort has been made to understand the flavour structure of quarks and 
leptons using non-Abelian family symmetries such as A±, S4, 50(3), SU(3), and so forth 
(see, e.g. pQ), where the three families of SM fermions can be unified in representations of the 
family symmetry. For reviews with extensive lists of references, we refer the reader to [2]. One 
strength of such non-Abelian family symmetries is that they can provide an explanation of 
the observed close to tri-bimaximal mixing in the lepton sector. Most of these models are also 
formulated in a SUSY context, although this is often not worked out in any detail. In many 
works only the superpotential relevant for the vacuum alignment of the family symmetry 
breaking Higgs fields, the "flavons", and the fermion masses and mixing angles is presented, 
without considering the consequences for the SUSY breaking sector. It is one main aim of 
this paper to highlight the entanglement between the SM and the SUSY sector when a careful 
model analysis is performed. 

Family symmetries (in particular non-Abelian families discussed towards explaining tri- 
bimaximal lepton mixing) can indeed constrain the structure of the sfermion mass matrices 
such that compatibility with FCNC constraints or LFV is improved, see, e.g. [3Tlll|. Several 
aspects regarding the connection between family symmetries and the SUSY flavour structure 
have been discussed in the literature: For example the role of induced flavon F-terms has 
been analysed in |12[ I13j . Canonical normalisation effects in family symmetry models in 
supergravity, and the dependence of their relevance on the messenger sector of the models, 
have been studied in [13]. In another recent study, the impact of the messenger sector on the 
predictivity of family symmetries for the soft SUSY breaking parameters has been highlighted 
[15J. Furthermore, in combination with family symmetries, supersymmetric Grand Unified 
Theories (GUTs) are an attractive framework for trying to solve the flavour and SUSY flavour 
puzzles, since they provide a link between the different (s)fermion species, i.e. between quarks 
and leptons respectively squarks and sleptons. 

Many studies in the literature focus either on explaining the flavour structure of the 
fermion masses and mixing parameters or on the SUSY flavour puzzle, but not on both at the 
same time. However, if SUSY will be discovered, it is clear that successful models of flavour 
will have to provide an explanation for both, the flavour structure of the Standard Model 
fermions and the flavour structure of their scalar superpartners. In this paper, we discuss 
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aspects of such "SUSY flavour" models, which aim at predicting both flavour structures 
simultaneously, in the context of supergravity (SUGRA). We point out the importance of 
carefully taking into account SUSY-specific effects, such as 1-loop SUSY threshold corrections 
and canonical normalization effects, when fitting the model to the data for fermion masses 
and mixings. These effects entangle the flavour model with the SUSY parameters and lead 
to interesting predictions for the sparticle spectrum as well as, for instance, for the neutrino 
parameters. In addition, family symmetries introduced to explain the flavour structure of the 
Standard Model fermions can also make predictions testable in future flavour experiments. We 
demonstrate these effects by analyzing an example class of flavour models in the framework 
of SU(5) Grand Unified Theories with a family symmetry with real triplet representations. 
Thereby we benefit from the GUT Yukawa relations y^/y s = 6 and y T /yb = 3/2 two of us 
proposed in |16j . From there we derive, in a semi-analytical as well as in a numerical analysis, 
the constraints on the SUSY parameters and the predictions for FCNCs, LFV and EDMs. In 
the class of SUSY flavour models considered as example, flavour violation through the SUSY 
soft breaking terms is dominated by the trilinear terms. We discuss the possible signatures 
of such a scheme of "Trilinear Dominance" El 

The paper is organized as follows: In Sec. [2] we first define the class of models we will 
consider as an example. At this point the model is not yet supersymmetric, i.e. we specify 
only the fields necessary to describe the SM fermion masses and mixing angles. In Sec. [3] we 
discuss the relevant SUSY specific effects, like canonical normalization and SUSY threshold 
corrections. There, we also discuss the structure of the soft SUSY breaking terms, first in 
the flavour basis and then also in the SCKM basis. The fourth section is devoted to a semi- 
analytical analysis, where we describe an approximate procedure to extract the low energy 
flavour observables from the high scale model. We also show how, from a fit of the low energy 
fermion masses and mixing angles to the data, constraints on the SUSY spectrum arise. In 
Sec.[5j we discuss the results of a MCMC analysis, which supports our semi-analytical results. 
We also present there phenomenological consequences for the SUSY spectrum and low-energy 
observables. The last section contains a summary and our conclusions. 

2 An Example Class of Flavour Models 

Before going to the SUSY specific ingredients we will first describe in general terms the class 
of flavour models we are going to consider as an example for all the SUSY specific effects. 
We emphasize that our analysis is not restricted to this special class of models. In contrary 
it should be applied to any given SUSY model of flavour to check its validity. As we will see 
it is necessary to include SUSY effects to check if a correct description of SM fermion masses 
and mixing angles can be obtained. 

2.1 Basic Structure of the Example Class of Models 

Our starting point for describing the class of flavour models, to be discussed as an example 
in the next sections, is a Grand Unified Theory (GUT) with gauge group SU(5) plus a non- 

5 We note that the scenario we consider here differs from the one proposed in |17) . where it is assumed that 
the fermion mixing is generated entirely by radiative effects especially from the trilinear terms. In our case, 
the SM and SUSY flavour structures are both generated from family symmetry breaking. Trilinear dominance 
means that the misalignment between the Yukawa matrices and the A-term matrices dominates the new physics 
flavour violation effects. 
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Abelian family symmetry group Gf with real triplet representations such as, for example, 
A4 or SO(3). For the matter fields of the SM fields S'C/(5)-unification implies that they are 
unified in the 5 representations Fi and in the 10 representations Tj, where % = 1,2,3 labels 
the three families, as follows: 
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We will assume that the three fields Fi are components of a triplet F under the family 
symmetry Gf- The other three fields Tj remain singlets under Gf- As we will see, this allows 
the top mass to emerge at renormalizable level, while the masses of the other lighter fermions 
will be generated from effective operators after family symmetry breaking. In addition to the 
SM fermions, we add two right-handed neutrinos Ni,N2 to the spectrum, which are singlets 
under the family symmetry, in order to allow for two massive light neutrinos via the seesaw 
mechanism as observed in neutrino oscillations. Note that this is a minimal setup and implies 
that one of the light neutrinos is massless by construction. It is straightforward to extend the 
model to include a third right-handed neutrino, to give mass to all three light neutrinos. 

As mentioned above, apart from the top quark mass all masses and thus also the mixing 
parameters of the fermions arise after the family symmetry (and the EW symmetry) gets 
broken. Following [18] we will assume that the family symmetry is broken by flavour-Higgs 
fields (the flavons) (pi in triplet representations of Gf, acquiring vevs in the following directions 
in flavour spac^J 



oc 
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The least understood and most difficult-to-test part of current flavour models with non- 
Abelian family symmetries is their so-called "vacuum alignment" sector. This is the part of 
the theory which addresses the question why the flavons point in these specific directions in 
flavour space. In this paper we will not enter this part of the flavour models in much detail. 
However, we will come back later to one (at least in principle) testable feature of the vacuum 
alignment sector in supersymmetric models, namely the size of the induced F-terms for the 
flavons. Finally, we remark that the scale of family symmetry breaking can in principle vary 
over a wide range. In the following, we will assume for definiteness that the family symmetry 
is broken close to the GUT scale. 

We will now describe the effective Lagrangian of the model class. For its origin from a 



renormalizable theory we refer the interested reader to Sec. 2.2 where we discuss the messenger 



field sector and sketch the additional discrete symmetries of the setup. In Sec. 2.3 we will 



6 We note that in explicit models one might probably split up (^2} into two flavons (cf. Appendix B of [18]). 
one with a purely real and another one with a purely imaginary vev. A mechanism for realising such vacuum 
alignments is presented in [T3] . 
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describe the resulting Yukawa matrices and the mass matrix of the right-handed neutrinos. 



A reader mainly interested in the phenomenological analysis may skip Sec. 2.2 



The effective Lagrangian of the class of models considered here has the following structure: 



-Cf =V2FH 2 J — ^-=-Ti + -^^T 2 + -^^T 3 ] H 5 + V2T 2 ^^FH 5 + H.c. , 



cff _./or ff J <h. rr_ , 02 ^ , 03 ^ fj , ,777^,^24 02 

M 2 M2 

(4) 

= J (itff + + ^2?) *5 + H.c. , (5) 

-£>t = f(^i + #5 + H.c. , (6) 



M 2 

-^. = ^ ? + ^NUM,.. (7) 



H 2 4 stands for a Higgs field in the 24 representation that breaks SU(5) to the SM, and 
H§ and stand for Higgs fields that contain the MSSM Higgses and H u in their SM- 
doublet components. The Higgs fields are singlets under the family symmetry. To shorten 
notation, 0(1) coefficients have been absorbed (without loss of generality) in the mass scales 
in the denominators, which correspond to the masses of heavy "messenger fields" as will 
be discussed below. As one can see, only ratios of flavon vevs over the messenger masses 
enter the Lagrangian, reflecting the fact that the family symmetry breaking scale is a priori 
undetermined. 



2.2 Messenger Fields and Discrete Symmetries 

Some properties of the effective Lagrangian are determined by the family symmetry Gf- For 
example, since F is a triplet under Gf it has to appear together with a flavon in the terms 
in Cy which lead to the down-type and charged lepton Yukawa matrices. This implies 
that the rows of Yd and the columns of Y e will be generated dynamically by the flavon vevs. 
Other properties of the effective Lagrangian, for example that the top quark Yukawa coupling 
appears at the renormalizable level while all other up- type Yukawa couplings do not, requires 
additional symmetries and the knowledge of the "messenger sector" , i.e. of the fields which 
are interchanged to generate the effective interactions. In summary, the class of models we 
like to outline here has the following features: 

• The rows of Yd and the columns of Y e will be generated dynamically by the flavon vevs 
with the rule that T\ couples only to <j>\, T 2 couples only to (p 2 and 4> 2 , and T3 couples 
only to </>3. 

• The Yukawa coupling ratios y T /yb and y^/y s at the GUT scale will be predicted as 3/2 
and 6, respectively, by Clebsch-Gordan coefficients from GUT symmetry breaking. As 
shown in [16] this requires a moderate or large value of tan f3. 

• At leading order Y u will be diagonal. 

• Left-handed neutrino masses will be generated via a see-saw type I mechanism with two 
right-handed neutrinos. 
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Figure 1: Supergraph diagrams inducing the effective superpotential operators for the down- 
type quarks and charged leptons. 



• The mixing in the neutrino sector is of almost tri-bimaximal form. 

• We consider spontaneous CP violation from the flavon vevs (see, e.g. [6]), in our 
specific case from (^2) only. 

The first point listed above is realized, for example, when the T, as well as the associated 
flavons carry some charge under a specific additional discrete symmetry, for instance a Z n 
symmetry. In addition, we have to make sure that the ratios among Yukawa couplings men- 
tioned above are reproduced by correct contractions of the SU(5) indices of the Higgs and 
matter representations. This can be achieved if the effective operators are generated by inte- 
grating out the appropriate messenger fields |16j . The relevant operators for the down-quark 
and charged lepton sector are, for example: 

^(^3)5^24(^5^)5 (8) 

leading to y T /yb = 3/2 and 

-l^(H^ 2 y 5 F(H 2 4T 2 )io (9) 

leading to y^/y s = 6 at the GUT scale. Here the subscript refers to the used SU(5) index 
contraction and M represents generically the mass scale of the messenger fields, which will 
be specified below. H24, and H5 are Higgs fields in the 24 and in the 5 representations 
of SU(5) respectively. H24 can be the usual Higgs responsible for the breaking of SU (5) 
to the SM gauge group by acquiring a vev (H24) = v 2 a diag(l, 1, 1, — 3/2, — 3/2). Clearly, 
the messenger fields, which generate these operators, can be distinguished by their charges 
under the additional discrete symmetries which distinguish between T 2 and T3. The diagrams 
generating the Yukawa couplings for down-type quarks and charged leptons are displayed in 
Fig. [TJ where the masses of the messenger fields are also shown. 

If we would introduce at least one discrete symmetry for each generation and the messenger 
fields generically carry charge under exactly one of these symmetries, it follows that the up- 
quark Yukawa matrix is generated with a diagonal structure. Specifically, we will arrange 
that only the top quark mass is generated at the renormalizable level, by associating T3 to a 
Z2 symmetry, for instance, while the charm quark and up quark masses are generated only 
after family symmetry breaking, by associating T 2 and T\ to larger symmetries as, e.g. Z n 
symmetries with n > 2. The relevant operator for the up-quark mass can then be 

^(TiTVjgtfsOMi)! , (10) 
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M [i) 




Figure 2: Supergraph diagrams inducing the effective superpotential operators for the up- 
type quarks. 




M (i) 





M, 



(i) 




Figure 3: Supergraph diagrams inducing the effective superpotential operators for the neu- 
trino sector. 



generated via diagrams as shown in Fig. [2} Similar considerations can be made for the neutrino 
sector, where now the RH neutrinos Ni carry the same charges as the Tj to associate them 
with the flavons. Majorana masses for the RH neutrinos and neutrino Yukawa couplings are 
then generated for instance by the diagrams shown in Fig. [3j 

In summary, we will consider models with auxiliary discrete symmetries to separate each 
of the Ti, leading to the following pairings between matter and navon fields: 



T± , N\ i — > (pi , 
T 2 ,N 2 ^ 4>i , 

T 3 <— > 03 • 



(11) 



We note that F is not paired to a particular navon since it is a triplet under the family 
symmetry and contains all three generations of the five-dimensional matter representations. 
Basically the following considerations are based on these pairings (and the sketched messenger 
sector) and do not explicitly depend on the discrete symmetries used to achieve them. In 
general these do not have to be chosen only to fulfill these pairings, but also to shape the 
navon potential such that it gives the desired vacuum alignment (cf . \18\ [19] ) . Since this task 
will require different quantum numbers for the flavons 4> 2 and it might be desirable to 
introduce another Higgs representation FL 2 t± or H^. Due to this difference we also assume that 
N 2 is not paired with (j> 2 . However, these model building issues are beyond the scope of this 
paper. 
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2.3 Resulting Yukawa Matrices and Right-handed Neutrino Masses 

The Yukawa matrices after GUT and family symmetry breaking are defined in the following 
way: 

-C Y = OGkQiUjH* + (Yj) lJ Q i d c J H d + (YZ)ijLie?H d 

+ (Y*)ijLii/jH u + (M v c) ijV fv c 3 + H.c. . (12) 

This convention matches the one used by the Particle Data Group |20| . In our class of models 
the Yukawa matrices and M v c have the following form at the GUT scale: 

'0 ei -ei \ 

Yd = | e 2 e 2 + i£2 e 2 + we 2 , (13) 
DO e 3 J 

ciei -ciei \ 

Y e T = ( C262 C2€2 + iC2e2 C262 + ^C2e2 J , (14) 

c 3 e 3 / 

<Vu 0\ 

K = | y c , (15) 

^/ 

y, 2 \ 

5^ = | Vv\ Vu2 , (16) 

-Vvl Vu2j 

— ( M o' ! (17) 

We remark that in the definition for the Yukawa matrices we have introduced a complex 
conjugation to match the PDG convention. This is the reason why +i instead of — i appears 
in the 2-2 elements of the down-type quark and charged lepton Yukawa matrices. 

The matrices Yd and Y e r are expressed in terms of the "expansion parameters" ex, €2, £2 
and £3, which are defined via the relations: 



V24 (01 ) / -, 1 v 24 (4>2) I 1 V24 




M[ Mi I / ' M 2 I / M3 M 3 



^24 



M 2 ' M 2 





(18) 



The structures of and , and its close relation, will lead to testable predictions of our 
class of models and the parameters ex, £2, e*2 and e 3 will be determined by a fit to the data 
for down- type quark and charged lepton masses and to the quark mixing parameters. The 
factors ci = C2 = c 3 = — 3/2 and c*2 = 6 are Clebsch-Gordan coefficients arising from GUT 
symmetry breaking and have been proposed recently in [16] as new possible GUT relations in 
gravity mediated SUSY breaking scenarios. Clearly, the obviously incorrect ££7(5) relation 
Yd = Yj would be obtained if all Clebsch-Gordan coefficients were one. To provide an idea 
of the size of the expansion parameters, let us say that in SUSY models with tan ft ~ 30 we 
will obtain small ei, e 2 and £2 of O(10~ 3 ) or less and e 3 of O(0. 1). 
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The parameters in the diagonal matrix Y u simply have a one-to-one correspondence to 
the up-type quark masses and the ratios of the squared elements of Y v over the right-handed 
neutrino masses in M u c stands, via the seesaw mechanism, in direct relation to the observed 
light neutrino mass splittings. By the structure of the vevs of <j)\ and $2 the leptonic mixing 
matrix is of tri-bimaximal form to leading order and receives calculable corrections by the 
(small) mixings of Y e , resulting in predictions for the neutrino mixing parameters as we will 
discuss later in the paper. 

Defining the ratios of flavon (and H24) vevs over messenger masses as 

^) = lv^- and C W = ^, (19) 

with V24 being the vev of #24, the entries of Y u , Y v and M v c are defined as 

M\ 2 , ^ M 2 2 2 „ 
Vu = MfMf^ 1 ' Vc = M^M' 2 U ' Vt = m ' ^ 

Vul = ^Cl > Vv2 = ^=$2 , (21) 

Now we have basically defined the building blocks of our class of flavour models. We 
would like to remark that within this class of models it is possible to define a specific version 
of a model by specifying the full set of (discrete) symmetries, the Higgs content at the GUT 
level, the vacuum alignment mechanism as well as the messenger sector of the theory. Note 
that, a priori, any family symmetry with real triplets can be used. Many predictions of the 
model class, however, are quite generic and will not depend on the realization in a specific 
model version. 



3 The SUSY Specific Parts 

In principle a model belonging to the class discussed in the previous section may also be 
realised in a non-supersymmetric setup. However, if SUSY will be discovered, it will be 
necessary to extend any given model of flavour to a supersymmetric one and to perform a 
careful analysis of its predictions in the SUSY sector. The first task, to be discussed in this 
section, is to specify the superpotential W and the Kahler potential K of the SUSY version 
of the model, which is consistent with the symmetries of the model and can be generated 
by its messenger field content. From W and K we can calculate the canonically normalized 
Yukawa couplings as well as the soft SUSY breaking sfermion mass matrices and the trilinear 
sfermion coupling matrices. 

In the following, we will focus on the scenario that SUSY is broken in a hidden sector and 
then mediated to the visible sector by gravity. We will also assume complete sequestering of 
K and W, meaning that we can split the superpotential and Kahler potential into two parts 

w = w vis + w hid , 

K = If vis + -Khid , 
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Figure 4: Supergraph diagrams inducing the effective Kahler potentiai operators for the 
fiveplet superfield F. 



respectively, where the visible superpotential is just W v i s = Wy d l + Wy u + Wy u + Wm v c + W 
and K v [ s = Kp + Kt + Kjj + where the dots refer to terms for the other visible sector 
superfields, and where W contains additional superpotential terms including, for instance, 
the flavon and GUT Higgs potentials. 

Furthermore, we will assume a gauge kinetic function such that unification of gaugino 
masses is realised at the GUT scale. 



3.1 The Superpotential 

The superpotential can be obtained by simply promoting the fields in the effective Lagrangian 
to chiral superfields. For our example class of models this give^] 

W Yit = V2FH 2i (j^Hl + -^T 2 + ^T 3 ) ft + ^f^ft , (23) 

wr '-i{*^ + jSv n+ ^ i3 ) B '- (24) 

w ^ =F {w Nl+ w N2 ) H "- (25) 

Promoting also the messenger fields to superfields, the effective operators in the superpotential 
can be generated by integrating out the heavy messenger superfields. Integrating out these 
messengers also induces effective operators in the Kahler potential, as we discuss now. 



3.2 The Kahler Potential 

Integrating out the messenger superfields generates not only the above superpotential oper- 
ators but also effective operators in the Kahler potential. An example, within our class of 
models, is shown in Fig. [4} The messenger superfield is the same that is interchanged in the 
left side of the left diagram in Fig. [TJ After family symmetry breaking, the kinetic terms 
are generically generated with an a priori non-canonical structure and have to be canonically 
normalized. This, in turn, has an effect on the fermion flavour structure, namely on the 



Yukawa couplings, as we will discuss in Sec. 3.3 



7 To simplify notation we will use the same letters for the fields also to denote the associated superfields. 
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From the diagrams that generate the superpotential terms, using the argument from above, 
one can easily infer the non-canonical Kahler potential operators which are generated by the 
messenger fields^ In our example class of models, the Kahler potential of the visible sector 
fields Tj and F, including these operators, is given by 

Kt =T lT }k<P + T 2 Tl ( kf> + ft(2)^4 N ) + T 3 Tlk'f , (27) 



K r =FFt (^>1 + 4 F) ^+4 fl ^+ if >^f) ■ (28) 

Note that Kp does not contain non-renormalizable operators with flavons, while such oper- 
ators are present in Kp and lead to non-diagonal terms after family symmetry breaking. In 
contrast, Kp does not contain non-renormalizable operators^] with -fi24- Gauge and family 
indices have been suppressed and are assumed to be contracted between the matter field and 
the flavon/GUT-Higgs field, see Fig. ^ The couplings l\ X ^ and depend on the fundamen- 
tal couplings of the messengers, but are essentially additional free 0(1) parameters as they 
do not appear somewhere else. 

From the diagrams in Fig. [T] and [2] one can see that there is also a flavon-dependent 
contribution to the Kahler potential of the Higgs field H%: 

KH-B&Jw+rpkjj). (29) 

We have not given here the Kahler potential for the right-handed neutrino fields since its 
effects are negligible for phenomenology or can be absorbed into other unknown quantities. 
We have also dropped terms induced by the couplings of the F and T fields to H§ and H^, 
since the resulting non-universal terms are suppressed by powers of the electroweak symmetry 
breaking (EWSB) vevs over the heavy messenger scale. 

3.3 Canonically Normalized Yukawa Matrices 

As typical in SUSY/SUGRA flavour models, the kinetic terms in our example class of models 
are generically non-canonical and have to be canonically normalized, which has an impact on 
the fermion flavour structure. We will now discuss the leading order effects of this procedure. 
The transformation to such a basis can be achieved as follows (see, e.g. [6]): 

F -> Pp l F , (30) 

with 

P F 1] K F P- 1 = 1 , or equivalently Kp = P F P F , (31) 

where Kp is the Kahler metric defined as Kp = 9 Kp + . Likewise, we define the transforma- 

dFidF] 

tion matrices Pp and Pg s which canonically normalize Kp and Kg . 



8 We note that, in addition, there may also be effects from Planck scale suppressed operators. With the 
messenger field masses sufficiently below the Planck scale, these effects are subdominant and can be neglected. 

9 This is only true up to dimension 6. There are dimension 8 operators involving both flavons and H24, 
which we, however, neglect. 
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As a consequence of the field redefinitions, the Yukawa matrices Y^ e are transformed as 
follows: 

Y d ^{P Q l ) T Y d P dc \ Y e (P^fY^ 1 . (32) 

The rescaling factor needed for canonically normalizing Kg 5 is readily absorbed into the 
definition of the e's. In contrast to Y d:£ , the Yukawa matrix Y u remains unchanged in form as 
one can absorb all of Pt in a redefinition of y u , y c , yt- 
Several remarks are in order here: 

• The parameters which are directly related to the fermion masses and mixings are the 
e's. On the other hand, in our example class of models the e's depend on the product of 
two ratios of flavon vevs over messenger masses, for example 63 ~ C(3)((3); an d similarly 
for the other e's as defined in the previous section. In the following, we will assume the 
case that and £(2) are comparatively large (> O(0. 1)). This in turn implies 
that, as we will confirm later by the fit to the data, that Q 2 )> C(i)j C 2 |j are very small 
such that we can neglect the corresponding canonical normalization effects in a leading 
order analysis. The only e parameter which will not turn out to be very small is €3, and 
therefore we will introduce both parameters ((3) and C 2 i explicitly. 

• Since Kt also involves couplings to the GUT breaking vev of H 2 ^ but Kp does not, we 
get the relations Pp = Pd c = Pl and Pq / P u c / P e c This modifies the CG factors 
between the three fields contained in the tenplet T 2 , i.e. C2 and c 2 , due to the vev of #24 
which treats those fields differently. However, since C 2 |, is very small, as stated above, 
this can be neglected. 

Under the assumptions specified above the main effect of canonical normalization stems 
from the change in Y d and Y e due to rescaling of the third component of the field F, which is 
parameterised by £3- With the definition (3 = (3(1^ /k^) 1 / 2 , and after absorbing all other 
rescalings in the e's, the Yukawa matrices have the following form at the GUT scale to leading 
order in ("3: 

(y u o\ 

Y u = y c , (33) 
\0 y t ) 

/0 ei -6i(l-KD \ 

Y d =\e 2 e 2 + ie 2 (e 2 + we 2 )(l - ^Cf ) , (34) 

\0 e 3 (l -Ka 2 ) / 

/ ciei - Cl ei(l- ^C 3 2 ) ^ \ 

= c 2 e 2 c 2 e2 + ic 2 e2 (c 2 e 2 + wc 2 e 2 )(l - 5C3) J . ( 35 ) 

V c 3 e 3 (l-K3 2 ) / 

/ y v2 \ 

Y v = y vX y u2 . (36) 

\-y vX (X-\ll) y v2 {\-\H)j 

3.4 Soft Breaking Terms in Supergravity 

Let us now turn to the soft SUSY breaking terms which can be calculated from a given K 
and W. Please keep in mind that we restrict ourselves to the case of hidden sector SUSY 
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breaking transmitted to the visible sector by gravity and that we consider for definiteness 
complete sequestering of K and W into hidden sector and visible sector parts. The hidden 
sector is assumed to contain chiral superfields which break SUSY by obtaining both a vev 
and a non- vanishing F-term of order Mpi and 777.3/2 Afpi, respectively. 

In supergravity before canonical normalization, the soft masses can be expressed in terms 
of the full visible sector Kahler metric K [21] : 

™lb = m l/2^ab ~ Frn{drndnKab)F n + F^(d m K ac )K~^ (d n K db )F n , (37) 

where n(m) runs over all (conjugated) superfields while b(a) only runs over matter superfields. 
Clearly, given the sequestering assumption, only visible sector superfields whose F-term is not 



vanishing contribute to the second and third term of Eq. ( 37 ) . In our class of models this is 
the case for the flavons and the GUT Higgs field, which acquire an F-term proportional to 
their vev [3]: 

F$ = C # m 3 / 2 ($> , (38) 

where the coefficients C^,., , Ch 24 could be computed after fully specifying the flavon and 
GUT Higgs superpotential. Given the "irreducible" contribution m 3 / 2 (<l>) to F$ [3], we can 
assume the C$ to be real 0(1) parameters^*} 

On the other hand the SUSY breaking trilinear terms [21] take the form: 

" "hi.: ,, f( R ^ 

|Whid| m \\ m M' pl 



A abc = e 2M pi j^F m i\d m -^\ d a d b d c w vis + d m d a d b d c w vis + ■■■ , (39) 



where subdominant terms have been omitted. However after SUSY breaking in the hid- 
den sector also the Yukawa couplings appearing in W v { s are no longer the physical Yukawa 
couplings for the low energy theory but take the form 



Miid 



w* 



Yabc = ^ , (40) 

|Whid| 



multiplied by the same factor also multiplying A abc . This means that Eq. (39) is simplified 
to 

A abc = F h y d h*^rj Yabc + F m d m Y abc H , (41) 

where h runs over hidden sector fields and m over visible sector fields. For this we used the 
fact that the rescaling factor only depends on hidden sector fields and thus commutes with d m 
due to sequestering. We can then define the usual mSUGRA parameter Aq = F^d^^f which 

depends on the unspecified hidden sector superpotential and Kahler potential. Regarding the 
Yukawa matrices we can absorb the factoipj exp ( ) W£ id / \ Wud \ into the previous definition 
of our Yukawa matrices. However, the expressions linking the parameters e's and £'s would 
be then modified by a 0(1) coefficient. 

From these formulae it is apparent why sequestering simplifies our construction for two 
reasons: 



10 The specific conditions for this have been examined in |12| and [13j . 

n ff we assume one hidden sector field X with Fx = 3mg/ 2 Mpi and canonical hidden Kahler potential we 
get as factor exp(^) = exp(^) (~ 2 if A = 2m ). 
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• It is known that models of the considered class with family symmetries with real triplets 
do not protect all particle species against large flavour violation in their soft mass matri- 
ces (see, e.g. [6]). In our case, while the five-dimensional representation is protected, the 
ten-dimensional remains unprotected. In particular, without sequestering there could 
be Planck scale suppressed couplings to hidden fields on the diagonal of the tenplet 
Kahler metric, inducing 0(1) splittings of the tenplet soft masses. This would generate, 
in the basis where is diagonal (the SCKM basis), 1-2 entries in the soft mass matrices 
of the tenplet constituents which are only Cabibbo suppressed. 

• With a sequestered Kahler potential all flavour violation effects come from the F-terms 
of the family symmetry breaking Higgs fields, i.e. the flavons, and in principle also 
from the F-terms of the GUT-Higgs fields. This clearly reduces the number of variables 
entering the soft terms and makes the model more predictive. 



3.4.1 Soft Masses 



From the expressions in Eqs. (27), (28), and (29), we can now derive the structure of the 



soft-mass matrices at the GUT scale. 

After canonical normalization the first term of Eq. (37) gives rise to a universal contribu- 
tion to the soft masses. However, due to the coefficients C$, Fm(dmd n K)F n is not rotated 
away by the same transformation making K — »■ 1. Thus, such a term gives rise to non- 
universal soft masses already at the GUT scale. The third term has the same effect, but 
provides only contributions oc £ 4 , which we neglect. 

If we take mo = TJ3/2J m t ne canonical basis, the soft mass matrices for the matter fields 
then result up to second order in £3 (neglecting all the small C's): 

mj = mgl , (42) 

4 



m}. = m 2 Q l, (43) 



'0 (P 

,0 1, 



ml = m 2 l-(C <j)3 ) 2 (!m 2 ( | , (44) 



and for the leptons 




m 2 L = m 2 l-(C^) 2 ( 2 m 2 [ 0] , (45) 

m| = mgl . (46) 

Additionally the Higgs boson H§ which contains the light Higgs doublet gets a non- 
universal soft squared mass 

4 = (1 - (<^ 2 ) 2 4Vo • (47) 



If we assume ^ to be 0(0.1) then the non-universality amounts to only a few percent, 
which is smaller than the accuracy of the 1-loop RGEs for the soft terms we use in our 
analysis. 

It is interesting to note here that due to the assumed pairing between different fields the 
soft mass terms are actually only diagonal for the tenplets. The soft terms for the fiveplets 
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have off-diagonal terms which are small since they are generated by terms of order the small 
Q 2 . The large £3 does not generate off-diagonal terms. However, this is only the case in the 
flavour eigenbasis. Upon rotating to the basis where the Yukawa matrices are diagonal (the 
SCKM basis) these non-universal diagonal entries induce off-diagonal entries suppressed by 
mixing angles. 



3.4.2 Soft Trilinear Couplings 



The first term of Eq. (41) does not give rise to flavour mixing, since it is proportional to the 



Yukawa matrix and hence diagonal in the same basis. 

The term proportional to d m Y has the same flavour structure as the corresponding Yukawa 
matrix, but it is not necessarily proportional to it, due to the flavon F-term coefficients C^.... 
Thus, the A-term matrices A d and A e will not be diagonal in the SCKM basis. As we will 
see, this represents the dominant source of flavour violation in our class of models. We will 
refer to this scheme as "Trilinear Dominance". 

After canonical normalization, the trilinear terms before going to the SCKM basis look 
like: 



Ah 



A Y d + C H24 m Y d 
( 



+ rn 



a 



V 



i 2 e 2 




C,<; 



-C^ei(l - \Cl) 
CV 3 e 3 (l - 



(48) 



At 



A Y^ + C H24 m Yj' 



+ mo 



3 C ( j >1 ei 
e 2 + 6iC^e 2 



A u = A Y U + m 



A v = A Y V + m 




1(7^61(1 - |C|) 
-|C0 3 e 3 (l - TjCf) 



tti) 



C(f> 2 yu2 

C<t> 2 ijv2 

Mi - K3 2 )/ 



(49) 



(50) 



(51) 



Note that the contribution proportional to Ch 2 ± is J us t Yd e times a factor and thus does not 



give any additional flavour violating effects, but only a splitting between Aq^ u and A, 
Thus for simplification, we neglect these terms and assume 



04/ e- 



Ch 2 < 







(52) 



for the rest of our analysis. The 3-3 element of the corrections to A u vanishes since we 
assume that for the third generation the messengers are absent as the top Yukawa coupling 
is generated on tree-level. A v is of less phenomenological importance as the right-handed 
sneutrinos get integrated out not much below the GUT scale. 



3.5 SUSY Threshold Corrections 

So far we have discussed the impact of a flavour model on the structure of the soft SUSY 
breaking sector at very high energy scales. However, to fit a SUSY flavour model to the 
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available data, there is a certain class of supersymmetric corrections at the low scale which 
is mandatory to include, in particular in the large (or moderate) tan B region. These are the 
well known SUSY threshold corrections |22j. 

For simplicity, in the semi-analytic part we will use the formulae for the tan /3-enhanced 
SUSY threshold corrections in the EW unbroken phase, where effects of electroweak symmetry 
breaking are neglected (see, e.g. [23J). We note, however, that it is important to include the 
SUSY threshold corrections for all three families of down-type quarks and charged leptons. 
The full formulae including EW symmetry breaking effects are given, e.g. in [24J. 

One-loop matching between the SM quantities and the MSSM ones at the SUSY scale 
-^SUSY leads to the approximate relations 

y^ r = (l + 6,tan^) ye M f T M cos^, (53) 
^ = (l + eg tan/?) yd M f M cos/3, (54) 
yl M = (1 + (e g + e A ) tan/3) y b MSSM cos/3 , (55) 



for the Yukawa couplings and 



nSM _ 1 + e g tan /3 ^mssm 
i3 l + ( e<? + e A )tan/3 i3 

nSM _ /iMSSM 
'12 — p 12 ! 



(56) 



(57) 

rSM _ cMSSM , K o\ 
°CKM — °CKM : l°°J 

for the quark mixing parameters. The one-loop SUSY corrections are parameterised in terms 
of the new e parameters. We see that although the e are loop suppressed the corrections can 
be sizeable if tan B is large enough. For explicit formulae and typical numerical values for the 
loop corrections, see [23], where we find 

tan B 

e/tan/3 = (ef + ef ) tan/3 = -(0.072 -r 0.020)— - , (59) 

oU 

e q tan B = (ef + ef + ef ) tan B = (0.030 -f- 0.248) , (60) 

oU 

e A tan B = e y tan /3 = sgnyl 4 (0.025 ^ 0.140) . (61) 

oU 

We should note that - apart from neglecting effects of electroweak symmetry breaking and 
sfermion mixing - this simplified parameterisation assumes that the corrections due to su- 
persymmetric gauge interactions are universal for all three generations. In general this is not 
quite correct since the sfermion masses entering the formulae are in general not degenerate. 
Even in the CMSSM where the masses are set to be equal at the GUT scale at low energies 
the third generation sfermions are lighter due to the effect of large Yukawa couplings on the 
running. In the full numerical analysis of Sec. [5j this will be taken into account. 

We will also see that these corrections enable us to constrain the CMSSM parameter space 
if we match our class of models to the experimental data at low energies. 

3.6 Flavour and CP Violating Effects 

The first thing to note is that in our class of models the leading order structure of the soft 
SUSY breaking scalar masses at the GUT scale is close to the one of the CMSSM. 
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A first potential deviation from the CMSSM flavour structure is represented by the non- 
universality of the third generation RH down squarks and LH sleptons given by the terms 



proportional to {C^ Z ) 2 C^ in Eqs. (44 45). This could in principle determine a misalignment of 
the quark (lepton) and squark (slepton) mass matrices already at the GUT scale. In fact, in 
the basis where the corresponding Yukawa matrices are diagonal (the so-called SCKM basis), 
m 2 ~ and m~ will acquire off-diagonal entries in the 1-3 and 2-3 sectors, approximately given 
by the product of the sfermion mass splitting and the corresponding rotation to the SCKM 
basis. For instance, we have: 



(2 ^SCKM _ / T ^i™2 



where is one of the two unitary matrices diagonalizing Y^: 



Y d 



V*Y d V$ . 



(62) 



(63) 



However, we see from Eq. (34) that the RH 1-3 and 2-3 rotations in are much smaller than 
the LH ones in V£ (which accounts for the CKM angles). Indeed, they even vanish in the 
leading order approximation (V^)i3 ~ (Yd) 13/ (Yd) 33. As a consequence, the source of flavour 
violation given by Eq. ( 62 ) should result subdominant with respect to the effects we discuss 



below. In the leptonic sector, it will turn out from the diagonalization of Y e (with numbers 
for the e's as given in Tab. [l| that (V#)23 ~ 10~ 2 . This can provide a source of r-fi lepton 
flavour violation at the level we discuss below. As we will see, such an effect is, however, of 
little phenomenological relevance, since r-fj, transitions result hardly in the reach of future 
experiments. 

Most notably, at the GUT scale the most relevant source of flavour violation is therefore 



represented by the trilinear couplings, Ad and A e , Eqs. (48 49). We will refer to this scheme 



as "Trilinear Dominance". The trilinear coupling matrices have a flavour structure which is 
similar to that of the corresponding Yukawa matrices, but they are not in general aligned 
to the Yukawas in flavour space. Such misalignment, as previously observed in the context 
of SU(3) [U El H2] and A4 [9j H3] flavour models, is a consequence of the flavon F-term 
coefficients C^.j, and it will clearly determine that the A-term matrices are in general not 
diagonal in the SCKM basis. Therefore, barring the unnatural case of all C</> (i) acquiring a 
common value, the A-terms will be a source of potentially large flavour violating effects in 
our class of models 12 Effects of large trilinear terms have also been studied in |17| 125]. 

Moreover, the large A-term matrices will induce flavour violating off-diagonal entries in 
the scalar mass matrices as well, through the RG evolution of these parameters from the GUT 
scale to low energy. This latter effect can be understood by means of approximate leading-log 
solutions of the 1-loop MSSM RGEs. Working in the basis of diagonal down quark and lepton 



Yukawas and neglecting for the moment the contribution of Eq. ( 62 ) , we have for the squark 
sector: 



(m 



^).™-^4)<>(^j 



(64) 
(65) 



12 On the other hand, we remark that depending on the flavon Higgs potential, the C^ (i . can also be strongly 
suppressed |12| . A discovery of SUSY particles at the LHC, but non-observation of any flavour violation effects, 
may thus point to a vacuum alignment with suppressed Cj,,.,, e.g. of the type discussed in |12| . 



17 



where (wi^)^ SSM represents the usual CMSSM contribution, proportional to entries of the 
CKM matrix. Keep in mind, that in the basis we have chosen there are only off-diagonal 
elements generated in m^, within the CMSSM. Analogously for the sleptons, we have: 

('4Wi = -^(^4feln(^), (66) 

K)^--i(4A)«ln(^IL), (67) 

where the subdominant SUSY see-saw contribution to (m?)j^j has been omitted. From the 
expressions above, we can see that at low energy we have sources of flavour violation both in 
the LH and in the RH scalar sectors, as well as in the left-right mixing sector (the A-terms 
themselves). What is interesting is that all these sources of flavour violation are determined 
by the matrices A& and A e , i.e. they are controlled by the limited set of free parameters C^ (i) 
(while the other parameters, such as the e^, are determined by the fit of the fermion masses 
and mixings), and the effects in the squark and slepton sectors are linked by the underlying 
SU{5) GUT. 

In order to give an estimate of the above discussed effects, we provide the following 
approximate expression of Ad in the SCKM basis (the expression for Aj is analogous with 
the substitution t{ — > Ci€i): 

A SCKM^ Ao Y d + mQX 

((c^ + C^-iC^fY-m (i^|-^) ei (i(C 13 + C 2 3u;|)-C 2 3U,) ei (l-iC3 2 )\ 
C 22 e 2 -iC 22 | + iC* 22 e 2 (C 23 e 2 + wC 23 e 2 ) (1 - 5C3) 

V °(§) -Ki) / 

(68) 

where Yd is the diagonalized Yukawa matrix and Cy = C^ (i) — C^,.. . 

From the expression above, we can see that, as expected, all the flavour-violating entries 
would vanish if Yd and Ad were aligned (i.e. = C^, 2 = = C^). Let us also notice 

that A^ CKM and A^ M acquire complex entries both on the diagonal and off-diagonal. This 
implies both flavour conserving and flavour violating contributions to the EDMs. 

We can now estimate the size of the flavour violating mass insertions (Mis), as usual 
defined as: 

J 

(4)W = T- (X,Y = L,R) (69) 
m f 

where rh 2 ~ is the average sfermion squared mass and a£ y are 3x3 blocks of the low-energy 
sfermion mass matrix. The off-diagonal flavour violating entries of A£ y read: 

^^ = K,>' (a£ r )^=^(4)^ ( A rl = a lr)> ( 7 °) 

where Vd is the down- type Higgs vev, Vd = v cos /3 (we consider in fact / = £, d). In the SCKM 
basis, the parameters (^xy)^j P rov ide a measurement of the misalignment among fermions 
and sfermions, responsible for the flavour changing processes. 
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Making use of Eqs. ( 64p7 ) together with Eq. (68), it is possible to obtain good estimates 
of the mass insertion parameters (S^ Y )iji including the CP violating phases. The results 
can be then compared with the bounds on mass insertions provided in the literature [26- 
129] . Such estimates also give a useful indication of the channels where the largest deviations 
from the SM predictions are expected and can be used to cross-check the results of the full 
numerical analysis (which will include the resolution of the full 1-loop RGEs, instead of using 
the leading-log approximation) presented in Sec. 5.2 Some relevant examples will be provided 
in Sec. 14. 51 



4 Semi-analytical Analysis 

After having discussed the GUT scale structure of the Yukawa matrices and the soft SUSY 
breaking terms we will now start with the discussion of low energy predictions of the considered 
class of models, first in a semi-analytical approach and then we will also present numerical 
results in the next section. 

We will start our discussion with some analytical estimates for the masses and mixing 
parameters at the GUT scale in the quark and the lepton sector which already show that 
our class of models is in principle realistic. Then we will also discuss the impact of the 
SUSY threshold corrections. Especially, we will discuss how these corrections can be used to 
constrain the CMSSM parameter space in our class of models. 



4.1 The Quark Sector 

We start with the discussion of the quark sector. Since the up quark Yukawa matrix is diagonal 
we can directly determine the CKM mixing parameters from the down quark Yukawa matrix. 
We find the following approximate relations at the GUT scale: 



qCKM 
'12 



gCKM 
'13 



qCKM 
7 23 



ei 



e 2 - i e 2 
£3 

e 2 + w e 2 



£2 



e 2 - i e 2 



Vd 



ei£2 



e 2 - ie 2 



£3 



Vs = |C2 - ie 2 | , 
Vb = e 3 (l - ^Cl 



(71) 

(72) 
(73) 



Note that we have to include the next-to-leading order term for #f 2 KM , which is comparatively 
large. Analogous terms for the other angles vanish due to the texture zeros in the third row 
of Yd- We can also derive an expression for the CKM phase Sqkm hi terms of the e's 



<5cKM 



arctan 



£2 
e 2 



(74) 



where we have dropped unphysical sign ambiguities in the e parameters. Actually the only 
physical signs are the relative signs between e 2 and e 2 and the sign of w leaving four possible 
sign combinations 
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For simplicity we consider the solution where all e parameters are 



positive except e 2 . 



13 However, requiring the correct CKM phase requires w > 0, which is already implied in Eq. (74 1 



19 



4.2 The Lepton Sector 

We discuss now the mixing angles and mass eigenvalues in the lepton sector. For the charged 
lepton sector we find at leading order 



'12 



7 13 



C2 £2 



e 



23 



C2 €2 

o, 
o, 



ic 2 e 2 



Ve 
Vr 



ci ei c 2 e 2 



C 2 e 2 — 1 C2 62 
C2 £2 — iC2 £2| 



(75) 
(76) 
(77) 



with the aforementioned ratios c\ = C2 = C3 = —3/2 and C2 = 6. Using also our formulae 
from the quark sector we can give relations for the GUT scale Yukawa coupling ratios 



Vr 

Vb 

lh 
Vs 
)k_ 
lid 



3 

2 ' 
6^1 

4y M 



15 



16 

0.4, 



cos 2 5ckm ~ 5.6 , 



(78) 

(79) 
(80) 



where we especially want to note, that the ratio y^/y s is not exactly equal to six as we 
naively imposed before, but is instead lowered by almost 10% due to corrections from €2- 
Nevertheless, all the ratios are still in good agreement with current experimental data |16| . 
For the neutrino mass matrix we find 



/ 



m,2 
IT 



where 7712 



2 C >3 



3_,2 K 
2»vl M v c ! 



2^3 




m 3 



(0 



v° 







\? 2 

2^3 



1P2 
2^3 



c 2 

S3 



\ 



(81) 



and ?7i3 = Vv2Wl~- The physical neutrino masses are then rri2 



^2(1 — |Cf) an d m 3 = 07,3(1 ~~ jCD- The lightest neutrino remains massless (mi = 0) in this 
minimal setup since we only introduced two right-handed neutrinos. 

Using the approximate formulae of |14| [30} [3"T] , we can obtain estimates for the mixing 
angles of the MNS mixing matrix as 



sin< 2 NS 



sin 



oMNS 
7 23 



qMNS 
'13 



^MNS 



1 
1 
1 

71 



1 



1 



( 



1 + ,ci 



yT+16 tan 2 (5ckm 
3 

7r + arctan(4 tan 5ckm) 



0.06 



1 + 4CI 



-7T 



(82) 
(83) 
(84) 
(85) 



where corrections to the neutrino mixing parameters due to £3 have been obtained using the 
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approximate analytical formulae from [31] , and sub-dominant corrections to #^ NS and <5mns 
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due to 8f 2 have been omitted. The model satisfies the generalized lepton mixing sum rule 
relations of [30] (which compared to |32| also include third family canonical normalization 
and RGE effects), which implies that the ^-induced correction to #J^ NS is small due to 
small cos(5mns)- RGE corrections are subdominant as well (and have been omitted) since 
the neutrino mass spectrum is hierarchical and since the neutrino Yukawa couplings are much 
less than one [33, 34J. In particular, we see that our class of models always implies MNS 
mixings angles larger than the tribimaximal values. The prediction for #^, NS , for example, 
will be tested at future neutrino oscillation experiments |35j. 

We also note that neutrinoless double beta decay and direct determination of the neutrino 
mass in this class of models is beyond the reach of ongoing and planned experiments due to 
the small neutrino masses implied by the hierarchy and m-i = 0. This could nevertheless 
be different in a variant of the model with a "type II seesaw upgrade" p2J [36] lifting the 
restriction mi = 0. 



4.3 GUT Predictions and SUSY Threshold Corrections 

Having the GUT scale structure of our model defined and using the simplified formulae for the 



SUSY threshold corrections as described in Sec. 3.5 we now derive constraints on the SUSY 
threshold corrections which we will then translate into bounds on the CMSSM parameter 
space. 



4.3.1 GUT Scale Yukawa Coupling Ratios 

We start our discussion of GUT scale Yukawa coupling ratios with the b-r Yukawas: 
3 _ y T (M Gl]T ) 1 + (e A + e q ) tan/3 y T (M GVT ) 



Vb(M, 



GUT J 



1 + qtan^ yb(M GVT ) 



(86) 



uncorrected 



where we refer with "uncorrected" to the GUT scale values without SUSY threshold correc- 
tions included. For simplicity we apply them here at the GUT scale, which is a good first 
approximation, see [23]. The value on the left side is the prediction of the model. 

The uncorrected GUT scale value of this ratio can easily be obtained, e.g. using REAP 
|33j. This yields the numerical value 



y T (M GUT ) 



y b (M GVT ) 



1.29 ±0.05 -0.03 



uncorrected 



tan j3 
30 



(87) 



where the uncertainty is largely dominated by the bottom quark mass uncertainty and the 
dependence on tan /3 is a very good approximation between tan /3 = 10 and 50. Plugging this 
into Eq. (86) we find the relation: 



(e A + e q ~ e t ) tan/3 » 0.16 ± 0.05 + 0.03 

Similarly, we find for the second generation 

V^Mgvt) l + e 9 tan/3 y M (M GUT ) 



tan /3 
30 



5.6 



y s (M GVT ) l + e,tan/3 y s (M GUT ) 



(89) 



uncorrected 
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However, in this case the uncorrected GUT scale ratio has a much larger uncertainty due to 
the large strange quark mass uncertainty. The value and range obtained from REAP is then 



y s {M { 



GUT J 



uncorrected 



4 3+ 1 - 7 



which can be translated into the relation for the threshold corrections 

(e 9 -e,)tan/3«0.32+°;l . 
Finally, for the first generation, we find 

y e (M GUT ) 1 + e 9 tan/3 y e (M GVT ) 



0.4 



with 



yd{M GVT ) 1 + e z tan/3 y d (M GVT ) 
y e (M G uT 



uncorrected 



Vd{M { 



gut; 



uncorrected 



which gives us the relation 



0-38±g;?f , 



(e g -e,)tan/3«0.05lg;|| 



(90) 



(91) 



(92) 



(93) 



(94) 



This constraint is even weaker than the one for the second generation and therefore does not 
give us new insights. Both constraints are not particularly severe. Comparing the relations 



to the ranges given in Eqs. (59)-(61) we expect that all constraints can be satisfied in a 
CMSSM-like setup. 

From these ranges we can also derive a rough estimate for the optimal value of tan /3 
needed to satisfy Eq. (88). Using the central values of the intervals of e q , €a and e; this yields 



tan ft « 20 , 



(95) 



from Eq. (88). In the more accurate analysis performed in |16| . we see that a b-r GUT scale 
Yukawa coupling ratio of 3/2 (within the CMSSM) can only be obtained for tan/3 ~ 30 or 
higher. For definiteness, we will fix tan /3 to the value 



tan j3 = 30 , 



(96) 



for our analysis. 



4.3.2 GUT Scale Mixing Angles 

For the mixing angles, we find equations similar to (186]) , namely 



eg KM (M GUT ) = (1 + e A tan/3) ^ KM (Af GUT )| uncorrccted • (97) 



Due to our matrix textures, we can always choose w using Eq. (73) to fit #03 M . However, 



ei and €3 are predetermined from the 1-2 sector and the third generation Yukawa coupling. 
To determine the model parameters, we can thus use the most constraining leading order 
equations, namely for y e , y„, 9f^ M , #23 and y T given in the previous two subsections. For 
this we use the numerical "uncorrected" GUT scale values and the approximation of applying 
SUSY threshold corrections at the GUT scale as before. The results are given in Tab. [TJ 
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Quantity Value 



e\ in 


10- 


4 


5.56 qi 


e 2 in 


10" 


-4 


6.36 qi 


e 2 in 


10- 


3 


-2.16 q t 


e 3 in 


io- 


1 


1.54 q T 



10 



2.66 g e + 0.30 



Table 1: Numerical values for the GUT scale model parameters. The q parameters are defined 
as qi = 1 + ei tan/3, q T = 1 + ej tan f} — A(f and = 1 + |(f + £a tan /3. 

For the "uncorrected" value of #f 3 KM , i.e. the one obtained without including SUSY thresh- 
old corrections, we find 

^ 3 KM (M GUT )| uncorrected = (3.24 ± 0.15) x 10- 3 , (98) 

and from Tab. [T] we get the model prediction at the GUT scale: 



#if M (M GUT ) = ^ = 3.62 (1 - ) x 10-^ , (99) 



where the dependence on e/tan/3 has dropped out. Plugging this into Eq. (97), we find the 
relation 

e j4 tan/3 + ^Cf ~ 0.106 ±0.043 . (100) 

We can remove £3 from this equation using the two neutrino mixing angles #^ NS and 
6>M NS , i.e. Eqs. p2l and (ph. This gives 



extant = 0.076lg;o53 • (101) 
In addition, using the values in Tab. [TJ we derive 

<5 C km = 1-28 , (102) 
which is consistent with experiment within ler. 

4.4 Constraints on the CMSSM Parameter Space 

So far we have derived constraints on the SUSY threshold corrections from the demand 
to match our assumed GUT scale structure of the Yukawa matrices to the observed SM 
quantities. 

In the next step we will use them to constrain the CMSSM parameter space that deter- 
mines the SUSY threshold corrections. This implies that we neglect the small deviations from 
the CMSSM in the considered class of models. This is a good approximation for all flavon 
F-term parameters except C^ 3 , which can enter the spectrum in non- negligible amount. 
However, it does so by lowering the right-handed sbottom and left-handed stau masses and 
changing the GUT scale Af, and A T trilinear coupling parameters. The first two enter the 
threshold corrections only through loop functions, while the last two do not enter at all or only 
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via running. Thus, barring extreme situations such as, for instance, large C^ 3 , the threshold 
corrections should be quite insensitive to Cfe. 



To test which parts of the spectrum satisfy the equations (88) and (101), we implemented 
full one-loop RG evolution of the soft term parameters from the GUT scale to the SUSY 
scale. Instead of a fixed SUSY scale Mstjsy we use the more common definition 

Msvsy = y/m^m^ , (103) 

where m^-. are the two stop masses. This scale is calculated via an iterative algorithm that 
repeats running between the GUT and the SUSY scale until Msusy changes less than 1%. 
At the SUSY scale the threshold corrections are included as described in Sec. 13.51 

Using this procedure we did a parameter scan for a CMSSM spectrum defined by the 
parameters mo, Mi/ 2 and Ao. As ranges we took mo G [0; 4000] GeV, Mu 2 G [100; 1500] GeV, 
ao = Aq/ttiq G [1; 3]. The motivation for the range of Aq/ttiq is the fact that the contribution 
from gauginos to the running of trilinear terms always pulls them to smaller or even negative 



values of At tx e^. Thus in order to satisfy Eq. (100) and minimize deviation from tri- 
bimaximal neutrino mixing, we need to have a large positive value for Aq. However, it is well 
known that Aq cannot be set to arbitrarily high value, to avoid charge and colour breaking 
(this requirement typically gives |^4o I ^5 3mo) [37]. Furthermore, we chose /i > due to 
the phenomenological requirements of the anomalous magnetic moment of the muon, while 
tan f3 = 30 was fixed as mentioned previously. 

For our selection of points, we drop any that do not satisfy the LEP bounds on superpart- 
ner masses [20] or do not trigger successful EW symmetry breaking, and also the ones where 



the lightest stau is the LSP, i.e. where it is lighter than the lightest neutralino. Out of the 



remaining set of points we simply select the ones satisfying (88) and (101). These are shown 
in Fig. [5| Note that the recent results from ATLAS [38] and CMS [39] only give a constraint 
for the few points in the lower left corner of the allowed parameter space. These constraints 
are also valid for large Aq/tuq and rather independent of tan/3, as discussed, e.g. in [40] . 
We would like to make a few remarks on the preferred parameter space: 

• We find that, generically, a value of ao > 2 is preferred. 

• For lower values of ao the results depend only weakly on mo and we need a small M\i 2 . 

• Larger values of M\i 2 in general need larger values of ao and show a stronger correlation 
between mo and M±/ 2 . 

• There is a lower bound on mo of about 600 GeV leading to a rather heavy SUSY 
spectrum. However, due to the large ao the third generation can still be relatively light. 

Note that, although at first glance our spectrum looks very similar to the ones necessary 
to achieve b-r Yukawa coupling unification in the large tan f3 regime [JT] , this is indeed not the 
case. There, ao = —2 (note that the different sign is crucial here) and furthermore, although 
they also need Mij 2 <C mo as we do, the hierarchy there has to be even stronger (see for 
example [i2H4"4] ) . 

Another interesting information is the range of achievable values for e^tan/3. From the 
set of selected points we find that the largest possible value is 0.16 and that only 2% of the 
set has e^tan/3 > 0.10. 
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Figure 5: This plot shows the sets of points satisfying the constraints in Eqs. (88) and (101). 
We only show the points with the least deviation for every given tuq and M 1 / 2 in our scan. 
The size of the points is scaled with the deviation, such that larger points agree better and 
smaller points agree worse. 



4.5 Estimate of Flavour and CP Violating Effects 



In the following, we make use of the expressions presented in Sec. 3.6 in order to provide an 



estimate of the main flavour and CP violation observables that will be subsequently studied 



numerically in Sec. 5.2 



4.5.1 LFV Processes 

Let us start considering \i — )■ e transitions. Within SUSY scenarios, the channel which 
currently provides the most stringent constraint to the [i-e sector is \i — > ej. Using the 



formulae given in Eqs. (66 67, 68) with ~ mo, the leading contributions to the [i-e mass 
insertions result: 

I05ll)2i| - ^C^(C h +a )c 2 e 2 c 2 e 2 ln^^- ~ 5 x HT 6 C 22 {C~^ + a ) , (104) 
I04O21I - ^C^d^c^ln^ x 10- 4 Ci 3 , (105) 



l(<%J 2 i| * —C 12 ci€i ~ 6 x 10- b C 12 , (106) 

v d „ _ _ . n . . -, n -6 f l TeV 



|(#r)2i| ^ -r~C 22 c 2 e 2 ~ 6 x 10^ — — C 22 , (107) 
mo V m o / 

where clq = Ao/mo. The numerical estimates have been made by means of the values of the 
parameters as obtained from the semianalytical analysis presented in the previous subsection 
(see Tab. [I]) and the running of the ^4g CKM entries has been neglected. As mentioned above, 
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the only uncertainty in the full determination of the flavour violating effects relies on the 
flavon F-term coefficients C^ ( .., which we assume to be 0(1) [4j. A comparison of the values 



in Eqs. (104 107) with the bounds on the mass insertions in [28, 29J shows that the dominant 
contributions in the \i-e sector are provided by the LR mixing insertions (5^ L )2i and (<5lr)21 ; 
i.e. directly by the off-diagonal entries of ^4g CKM . This makes it particularly simple to estimate 
the resulting fj, — > cy rate. 

On general grounds, the branching ratio of ii — > ij'j can be written as: 

BR(£ t ^£ j7 ) _48vr 3 a n/lL|2 , . . R|2 



^r + l^r), (108) 



BR(k -> ijViVj) " G 2 F Vl l]l 1 13 
where the amplitudes A L , ^4 R , in terms of LR mass insertions only, approximately read 

MM- 

4^S^~ 2 ^ 4 = 4 ( L ^ R )- ( 109 ) 

J 4-7T m~ J J 

Here is the ith generation lepton mass, M% is the Bino mass, x\= Mf/m- and the loop 
function f(x) reads: 

v l-5x 2 + 4x + 2x(x + 2)lnx 

/(*) = • 



Inserting the values of (5^)21 and (^lr)21 provided in Eqs. (106, 107) in the above formulae, 



we can finally get the following expression (which neglects subdominant contributions): 



lTeV\ 4 / Mi 



BR(^e 7 ) = 3.5X10-^+^— j [^f) , (111) 

where we assumed M\ <C m; ~ rriQ, a condition which is often verified within the model, as 
we will see in the next sections. Remarkably, BR(^ — > ej) can be easily in the reach of the 
MEG experiment at PSI [36], whose final sensitivity is 0(1O -13 ). On the other hand, we see 
that just a conspiration of the CV 4 , coefficients can result in a value of BR(^ — > ej) at the 
level of the present experimental limit (1.2 x 10~ n )p*| 

Another promising observable in the /i-e sector is represented by fi — > e conversion in 
nuclei. Within SUSY models, the latter process is typically dominated by the same penguin 
diagram contributing to /i — >• &f (while a possible box diagram contribution turns out to be 
negligible) with the photon attached to a proton line. As a consequence, \i — >• e conversion in 



nuclei just requires an additional electromagnetic vertex. We can then make use of Eq. ( 111 ) 
to easily estimate the conversion rate: 

CR(/i -> e in N) ~ a cm x BR(p -»• erf) . (112) 

Let's now consider the LFV in the t-\jl sector, namely r — > fij. The approximate expres- 



In the numerical analysis presented in Sec. 



5.2 



to ICjjjl, | C 2 2 1 < 6. Thus, according to Eq. (Ill I, we can get at most BR(^i — > &y) ~ 2.5 x 10 



we varied |C^,,„ | in the range 0.3 -f- 3, which corresponds 



-ii 
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sions for the T-fi mass insertions read: 

\(SL)32\ * (^ 3 C 3 ) 2 ra 2 3 + ^C^ 2 +a )w(c 2 ~e 2 ) 2 In 

~ 10- 2 (C </)3 C3) 2 + 2.5x 10- 4 C2 3 (C^+a ), (113) 
l( 5 RR)32l - ^23(^3 + an)w c 3 e 3 c 2 e 2 In ~ 8 x 1(T 3 + a ) , (114) 

I^R L )32| - ^-Cfc wc 2 ~e 2 ~ 2.5 x 10- 4 (l^f) q 3 , (115) 
mo V m o / 

K^r) 32 | * o(— WOaal * 0.051(5^)32] - (116) 

From these equations, we see that the main contribution to |(<5f L )32| ~ {C^^) 2 ^^)^ does 
not arise from the running, but directly from the non-universality of the third generation 
LH sleptons, as discussed at the beginning of Sec. |3.6| All the mass insertions can give a 
contribution to the r — > fij rate of a comparable size, so that the situation is more involved 
than in the fj,-e sector. The total rate of r — > ^7 should, however, be too small for detection at 
future experiments. For example, considering only (6^)32, we can write the A\\ amplitude 
like the following |47| : 

^ oil^M (n7) 
32 4vr 60 ml V ' 



Using this expression with the mass insertion as in Eq. (114), Eq. (108) then gives: 



BR(r M7 ) ~ 1.6 x lO^Cf, (C 03 + a ) 2 (^) , (118) 

which gives at most BR(r -> ^7) ~ ©(lO" 10 ). The other mass insertions give slightly 



smaller but comparable contributions. As we will see in Sec. 5.2, a full numerical evaluation 
shows that the different contributions can sum up to somewhat larger values of BR(r — > fij) 
(< 0(1O~ 9 )), but still barely within the sensitivity of future experiments. 

In the r-e sector, the rates are even smaller. For instance, we find BR(r — > e"f) < 10~ n , 
far beyond the reach of any foreseeable experiment. 

4.5.2 Electric Dipole Moments 

In a similar way we can estimate the size of the electric dipole moments as predicted by 



the model. From Eq. (68), we see that the A-term matrices, and A e , acquire a CP 
violating phase in the 1-1 entry, again as a consequence of the misalignment with respect to 
the corresponding Yukawas, induced by the flavon F-term coefficients. It is well known that 
complex A-terms can easily induce unacceptably large EDMs (this is often regarded as part 
of the "SUSY CP problem"). It is therefore remarkable that, in our class of models, not only 
the phases of the flavour changing entries but also the flavour conserving ones get suppressed 
by the expansion parameters e«. For instance, we see that the imaginary part of (yl^ CKM )n 
acquires an additional e 2 /e 2 suppression (i.e. of the order of the Cabibbo angle) with respect 
to the real part. This is even more effective in the leptonic sector, where we have (c 2 e 2 ) / (c 2 e 2 ) , 
with c 2 = —3/2, c 2 = 6. Complex A-terms clearly induce complex mass insertions through 



the RG running (see Eqs. (64)-(67)). Therefore flavour mixing induced sources of EDMs will 
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be also present within our model 15 Even though we are going to take into account all these 
contributions in the numerical analysis, let us here provide an estimate of the electron EDM, 
as induced by a complex (#lr)h only. 

In the regime Mi <C fhf, we can approximate the electron EDM like the following 



aiMi, 
An fh?~ 



11 



1119) 



where the imaginary part of (#lr)h can °e obtained from Eq. (68): 



LRjll 



— —C. 



22 



(C 2 £2j Cl£l 
(C 2 6 2 ) 2 



(120) 



Therefore, we get: 



- ~ 2.2 x 10" 29 



Mi 



100 GeV 



lTeV\ 3 

Coo cm . 

m J 22 



fl21) 



which gives at most d e ~ 0(1O -28 ) e cm, one order of magnitude below the present exper- 
imental limit. Therefore, the model predicts the electron EDM at a level that is at present 
consistent with non-observation, but still several orders of magnitude above the SM predic- 
tion, such that it can be easily tested by the currently running experiments. 
A similar situation is verified in the hadronic sector, where we have: 



(5 



LRjH 



Vd_ c ejei 

,yy~i _ 22 ~2 



3 x 10 



-7 



1 TeV 



C. 



22 ' 



:i22) 



which is about one order of magnitude below the present bound 
atomic EDM searches. 



, given by neutron and 



4.5.3 CP Violation in the Kaon Mixing 



Let us finally consider the bounds coming from SUSY contributions to the Kaon mixing, which 
are the most stringent ones within the hadronic sector. As usual, the strongest constraint 
is given by the CP violation parameter ex- In our model, all the relevant phases can be 
derived from the expression for ^ CKM given in Eq. (|68h . As in the leptonic sector, the 



largest contributions come directly from the off-diagonal entries of the A-term matrix. The 
main contributions to ex are then given by: 



(#LR.)l2 

(<5rl)i2 



Vd 


^22 


£2£l 








fflq 




Vd 


^22 


4 


rhg 


£2 



~ 10" 



1 TeV\ 

) °22 ' 



ma 



1 TeV^J 



rrin 



^22 ■ 



(123) 
(124) 



The resulting values are between one and two orders of magnitude below the present exper- 
imental bound, ~ 10" 4 [27]. Still, as we will see in Sec. 5.2l the interplay among different 



contributions can provide sizeable deviations from the SM prediction for ex- 



For a recent discussion of flavour violating contributions to EDMs, see for instance |48| . For studies about 
flavoured CP violation within 517(3) flavour models, see [8]. 
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5 Numerical Analysis 



5.1 Markov Chain Monte Carlo Analysis 

In the previous section we have discussed simple semi-analytical estimates, which showed that 
our class of models first of all seems to be compatible to current experimental data for fermion 
masses and mixings and second we have demonstrated how the CMSSM parameter space can 
be constrained. 

In this section we will improve our analysis using a Markov Chain Monte Carlo (MCMC) 
analysis. From this we can derive posterior probability distributions for the model parameters 
based on the experimentally measured quantities and uncertainties. Furthermore, it allows 
us to incorporate the deviations from the CMSSM in the considered class of models and fit 
the resulting spectrum to the data. Finally we will discuss predictions for various observables 
derived from that. 



5.1.1 Procedure 

The MCMC analysis was done using a standard Metropolis-Hastings algorithm [20] . For the 
conditional probability density in the usual form P(x\9) 

F(^)ocexpf-^V (125) 



where x represents measurements and 9 model parameters, and where with x 2 we denote the 
usual sum of deviations squared for the fitted parameters normalized with the experimental 
uncertainties. The quantities we use are the down-type quark and charged lepton Yukawa 
couplings, CKM mixing parameters and the 1-2 and 2-3 MNS mixing angles. Consequently, 
we varied only the ei parameters, C3, w, C^ 3 and the SUSY breaking parameters mo, M1/2 and 
ao = Aq/tuq. The other parameters, i.e. the up-type quark Yukawa couplings and neutrino 
parameters were fixed to their respective GUT scale value obtained neglecting SUSY threshold 
corrections. 

For the experimental values the quark and charged lepton masses at mt{mt) were taken 
from |49| as input. The quark mixing angles were taken from the global fit parameters as 
given in [20]. For the parameters in the neutrino sector we used the values given in [50J. In 
our analysis, we take the charged lepton Yukawa coupling uncertainties to be 1% (although 
they are of course measured with higher precision). A higher accuracy would not be justified 
as we are using only one-loop RGEs, for instance. 

These low energy observables are matched to the underlying model parameters to obtain 
the desired x 2 y i a the following procedure: 



1. Take the Yukawa matrices at the scale Mqut as shown in Eqs. (33)- (35). 



2. Evolve the Yukawa matrices from Mgut to Mstjsy using the one-loop MSSM RGEs. 
Right-handed neutrinos are integrated out at their mass scale with the appropriate 
matching as given by the see-saw formula. 

3. Match MSSM and SM Yukawa couplings incorporating the SUSY threshold corrections 
as described in the previous section. 

4. Evolve the Yukawa matrices to the top mass scale mt{mt) using SM RGEs. 
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Quantity 


Value 


Number of Chains 


5 


Number of Distinct Points 


468238 


Number of Iterations 


7435125 


Minimal total \ 2 with n^of = 2 


2.92 


PSRF - 1 


4 x 10" 4 



Table 2: Statistical characteristics of the Markov chains. PSRF is the maximal potential 
scale reduction factor [51] over all parameters. 



For the prior probability densities we used flat distributions for almost all varied param- 
eters. Of these the following are notable: The massive parameters mo and were taken 
from the interval [0; 4000] GeV and [150; 1500] GeV respectively. The prior for the trilinear 
coupling parameter ao = Aq/ttiq was taken be flat on the interval [0; 3]. Furthermore, we took 
Q instead of (3 to be flatly distributed due to the observation that the relevant SM quantities 
depend linearly on the square of £3. Thus a Gaussian distribution for these quantities also 
implies one for Q (neglecting boundary effects). The flavon F-term parameter C^ 3 was taken 
to be flatly distributed on the interval [—3; 3]. 

Based on this we calculated several Markov chains with the characteristics shown in Tab. [21 



5.1.2 Results 

The standard interpretation of the points found by the Markov chains is that they are a 
sample from the posterior probability distribution P(9\x). We can thus give probabilities for 
the model parameters and therefore 68% intervals and multidimensional boxes, where the 
parameters lie with that probability. The values obtained from our analysis are shown in 
Tab. [3| We did not show values for the C^ 3 and massive parameters tjiq and M±/ 2 as they all 



have a rather flat distribution, i.e. excess kurtosit 16 < —1, over the whole range allowed by 
the prior. However, the Markov chains prefer heavy masses and disfavour large positive Cs 3 . 

Based on this sample, we can also derive a probability distribution for each fitted SM 
quantity, interpreting them as functions and thus new random variables dependent on the 
model parameters. In the best case scenario they should be centered exactly on the experi- 
mental input value with a standard deviation matching the experimental uncertainty (or the 
one put into the analysis) and naturally be normally distributed. Deviations from this can 
naturally arise through model inherent correlations, prior effects and lacking convergence of 
the MCMC data distribution. The quantities obtained from our analysis are shown in Tab.|4| 
Looking at those numbers, we see that the considered class of models can in general 
describe the measured SM quantities quite good. The largest deviations are a very pronounced 
value for Sqkm at the la interval boundary, which can be attributed to the highly constrained 
1-2 sector of the Yukawa matrices and Y e . Then there is a shift of the value of #^ 3 KM also by 
about lcr and the neutrino mixing angles by la, each however with large standard deviations. 
This is expected from the range of possible values for ca together with the form of the 
constraint given by Eq. (100), which forces us to choose between a deviation in 0f 3 KM and a 



deviation in the neutrino mixing angles. Analogously, this also explains the lower value for 



3 For comparison, we note that a completely flat distribution would correspond to an excess kurtosis of —1.2. 
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Parameter 


Value with 68% Range 


ei in 10" 4 


5.88 ±0.11 


e 2 in 10~ 4 


6.64 ±0.13 


e 2 in 10~ 3 


2.27 ±0.04 


e 3 in 


1.673 ±0.073 


It 2 

2 S3 


0.035 ± 0.033 


W 


3.14 ±0.12 


ao 


2.13 ±0.51 



Table 3: Distribution characteristics of model parameters as found by the full numerical 
MCMC analysis. 



Parameter Experiment MCMC 68% Range 



0CKM 

oCKM 

"23 

0CKM 

<5cKM 




0.2277 ±0.0010 
0.0414 ±0.0012 
(3.60 ±0.17) x 10~ 3 
1 on9+ - 078 

i -^ uz -0.043 


0.2276 ±0.0010 
0.0414 ±0.0012 
(3.76 ±0.12) x 10~ 3 
1.286 ±0.003 


Vd in 10" 
y s in 10" 
Vb in 10" 


-5 
-4 
-2 


1.59 ±0.68 
2.99 ±0.86 
1.579 ±0.051 


1.62 ±0.53 

2.63 ±0.86 
1.564 ±0.047 


He in 10" 
Vn in 10" 
y T in 10" 


-6 
-4 

-2 


2.78 
5.88 
1.00 


2.78 ± 1% 
5.84 ± 1% 
1.00 ±1% 


flMNS 

ziMNS 

"23 

^MNS 

<>MNS 




0.602 ±0.017 
0.749 ± 0.071 


0.625 ±0.008 
0.825 ±0.018 
0.053 ±0.001 
4.81 ±0.05 



Table 4: SM quantities obtained from the MCMC data. The minimal \ 2 (including down- 
type quarks and charged lepton Yukawa couplings, CKM parameters, 6>^j, NS and 6>^ NS ) is 2.92 
with two degrees of freedom corresponding to a p- value of 23%. 
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Figure 6: Points found by the Markov chains with x 2 < 5.99 (p-values larger than 0.05). 
Points are ordered with respect to their x 2 such that the points which are fitting the experi- 
mental values better are drawn on top. In addition, point size scales with % 2 , making points 
agreeing better with experimental data larger. 



w compared to the value given in Tab. [Tl since it does not have to counter- act such a large 

~ 2 — 

e^tan/3 + 5C3 to arrive at an acceptable #23 KM . 

Another possible interpretation of the Markov chain data is that it just constitutes a 
sophisticated scan over the parameters that happens to be more fine-grained in the interesting 
regions of low x 2 ■ We can thus select from the chains those points that yield a x 2 lower than 
some threshold value and demote the Bayesian probability interpretation to the nice side-effect 
of enhanced resolution in parts of the parameter space 17 



5.1.3 Constraints on mo, M1/2 an( l Aq 

Based on the full MCMC analysis, we can also compute the constraints on mo, Af 1/ / 2 and A$, 
in analogy to the plot shown in Fig. [5] where we used the analytical approximations. The set 
of points with x 2 < 5.99 is shown in Fig. [6j With the two residual degrees of freedom of our 
fit this corresponds to p-values larger than 0.05. 

We can see that both spectrum plots look fairly similar. One difference, however, is that 
in Fig. [6] we find light orange as well as light blue points in the central region of the points, 
whereas in the plot in Fig.[5]there are mainly light orange points. This is explained by the fact 
that the colour of the points in Fig. [^indicates the value of Ao/mo which provides the best fit 
- and Aq/itlq = 2.25 turned out to be slightly better in this respect than Aq/tjiq = 2 - whereas 
in Fig. [g] all points of the MCMC sample with x 2 < 5.99 are shown. Another difference can 
be seen in the lower bound on rriQ. It moves from 600 GeV to 1 TeV due to the fact that 

17 Strictly speaking this is not necessarily true because the MCMC algorithm can stay at the same point 
for more than one iteration resulting in a possibly crucial difference between point and probability density 
distribution. 
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the points on the edge of the allowed parameter space of Fig. [5] need larger values of 1 /2<f| 
which are allowed by the approach used for the scan but disfavoured by the experimental 
values of the neutrino angles. This leads to an exclusion of the points with lower mo by the 
requirement x 2 < 5.99. 

Thus we can conclude that our approach in the simplified parameter scan is useful to get 
a first impression of the allowed parameter space. However, we stress that a full numerical 
MCMC analysis is still necessary to obtain accurate constraints, e.g. for mo- Furthermore it 
carries additional useful probabilistic information, e.g. in the apparent point density. 



5.2 Phenomenological Consequences 

In this section, we present the results of a numerical analysis performed by scanning the 
parameter space favoured by the MCMC fit presented above. The SUSY parameters mo, 
M1/2, A)> a s well as eu\, w, £3 and were taken from the MCMC fit discussed in the 
previous section. We remind that for the fit we used tan f3 = 30. The coefficients relating 
flavon F-terms and flavon vevs were then randomly varied in the following ranges: 



0.3 < IC^I < 3, 0.3 < IC^I < 3, 0.3<|C^|<3. (126) 
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Following the discussion of Sec. 3.3, the parameters (fi have been assumed to be 0(1), 



with the exception of £| 4 which has been set to 0.1. 

In order to compute the SUSY spectrum, we numerically solve the full 1-loop RGEs of 
the MSSM (with 2-loop RGEs for gauge couplings and gaugino masses) from Mqut down to 
the SUSY scale -Msusy = ■sj™l x m t 2 \\ For each point of the parameter space, we impose the 
following requirements: (i) successful EWSB and absence of tachyonic particles; (ii) limits on 
SUSY masses from direct searches; (iii) neutral LSP. Then, we compute the leptonic processes 
by means of a full calculation in the mass eigenstate basis |52j, the hadronic processes by 
means of the mass-insertion approximation formulae in [261 153|. and the BR(i? — > X s ^f) using 
SusyBSG [53]. We then require that the resulting BR(£? — > X s ~f) does not deviate from the 
experimental value [55J in more than 3cr (2.75 x 10~ 4 < BR(5 X s j) < 4.25 x 10~ 4 ) and 
thatBR(5 s fi + fT) < 5.2xl0~ 8 , 0.77 < \e K /e c ^ p \ < 1.23 [26l[56]. The LEP limit m^ > 114 
GeV does not impose additional constraints to the parameter space. 



5.2.1 Spectrum 

The MCMC sample we used for this analysis corresponds to the ranges for the SUSY particle 
masses reported in Tab. [5} The lower values shown in the table give an indication of the 
SUSY spectrum which can better account for the fermion masses and mixing: in the case of 
sfermions they correspond to the fact that for a given value of -Mi/ 2 a minimum value for 
mo is selected by the fit (as depicted in Fig. [6]). On the contrary, the upper values are just 
related to our choice of scanning the parameter space up to mo = 4 TeV, My 2 = 1.5 TeV. 

As we can see a rather heavy spectrum is selected, e.g. the heavier squarks are always 
above the TeV scale. However, gluinos and lightest stop and sbottom (t\, b±) can be much 
lighter. In particular, is driven to light values by the large value Aq > 2mo required by 

18 We notice that a full 2-loop computation might be necessary for a correct evaluation of the spectrum, 
especially for regions of the parameter space where EWSB requires precise tuning. However, we checked that 
for most of the parameter space 1-loop RGEs with the inclusion of 2-loop contributions to the ft functions of 
the gauge couplings and gaugino masses is accurate enough for the estimates we present here. 
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m s (GeV) 

Figure 7: Squark (blue: gx; green: b\; red: t\) versus gluino masses from the MCMC analysis. 
See the text for details. 



Particle Mass Range (GeV) Particle Mass Range (GeV) 
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107 4- 


1267 


552 4- 
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383 4- 


3680 


1077 4 


-4807 


£ 


1048 4 


-4091 


115 4 


- 129 


A 


261 4- 


4480 



Table 5: SUSY mass ranges used in the analysis. The upper values just reflect the choice of 
scanning up to mo = 4 TeV, M1/2 — 1.5 TeV. On the contrary, the lower values result from 
the successful fit of fermion masses and mixings and/or the applied constraints (see the text 
for details). 



the fit. In Fig. [7J we show the heavy squarks mass (blue) and the third generation squark 
masses (green) and (red) versus the gluino mass m g for the scanning of the parameter 
space described above. We see that the squarks turn out to be always heavier (even much 
heavier) than gluinos in the SUSY parameter space region favoured by our model, with the 
possible exception of ti, especially for large values of M\m (i.e. of nig). 

Even if gluinos might be light enough for an early SUSY discovery 19 it is clear that such 
an eventuality would select a corner of the parameter space (the bottom- left one in Fig. [6]). 
It is remarkable that the spectrum itself gives the possibility of testing our class of flavour 
models and distinguishing them from other scenarios. In this sense, the heavy spectrum, the 
squarks heavier than gluinos and the relatively light third generation are all generic features 
of our scenario, which the LHC could test after enough years of data taking. 

Sleptons are usually heavy with the possible exception of the lightest stau. In fact, we 
recall that most of the parameter space points lie in the region where mo 3> M\/2 and 
thus mi ~ rrig, since relatively light gluinos do not drive the squark masses to much larger 

19 With 1 fb _1 of collected data, the discovery potential of the 7 TeV LHC run, might reach : 
in the case, like in our scenario, of m§ <C rriq [57| . 
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values than sleptons in the RG running. In particular, the slepton masses turn out to be 
always larger than the mass of the W (and thus of the lightest chargino and the second 
lightest neutralino). This means that they cannot be efficiently produced in cascade decays 
of squarks to neutralinos. On the other hand, the direct Drell-Yan production cross-section 
for rrig > 200 -j- 300 GeV might be too small to have a slepton signal above the SM background 
[58]. As a consequence, sleptons seem to be hardly detectable at the LHC, within our scenario. 
However, the lightest stau represents an exception also in this case, since in some regions of 
the parameter space < m^o and the decay X2 ~~ ^ ~~ y Xi TT is open and potentially 
observable at the LHC. 

In our scenario, not only a specific region of the CMSSM parameter space is selected, 
but also effects from may result in non-negligible differences compared to a CMSSM 
spectrum. Observation of such a C^ 3 signature would of course be interesting, however, since 



only the sbottom and the stau masses are affected (cf. Eqs. (44 45) and Eqs. (48 49)), this 
seems to be experimentally challenging. 

The LSP is mostly Bino over all the parameter space. This is because the large At re- 
quired by the fit drives the running of m 2 Hu to large negative values so that the Higgsino 
mass parameter /i results always quite large (in particular [i 3> Mi), once the EWSB con- 
dition, fi 2 ~ —rrijj, is imposed. In fact, we find fi > 650 GeV over all the parameter space 
and consequently a quite small Higgsino component of the lightest neutralino (< 0.1). As a 
consequence, the LSP annihilation cross-section is not enhanced by a sizeable Higgsino com- 
ponent, and the so-called focus-point region for neutralino dark matter [59j is not viable in 
this scenario (for tan/3 = 30). On the other hand, there are some points of the parameter 
space (giving — m ^o)^ wnere an efficient stau- neutralino coannihilation [60] can reduce 
the neutralino relic density within the WMAP bound. This is because, as noticed above, 
the lightest stau fi can be driven light for sizeable values of Aq and C^ 3 , for which A T gets 
increased (cf. Eq. M9l). Moreover, the non-universal contribution to ^jy d , shown in Eq. (47), 
can be used to lower the mass of the CP-odd Higgs, A, triggering the so-called A-funnel region 
(corresponding to ttla — 2m^o), where the LSP annihilates efficiently through an s-channel 
A exchange [61] ■ In our model, this opens up the possibility to access the funnel region even 
for moderate values of tan/3. 

5.2.2 Flavour observables 

Let us now consider the predicted rates for the LFV observables fj, — > &y and r — > ^7. In 
Fig. [8j we plot BR(/x — > ej) versus BR(t — > ^7), for the variation of the parameters (in 
particular the C^... ) described above. We see that our class of models is at present not much 
constrained by the experimental limit on BR(/x — > cy), which is shown in the figure as a black 
thick line. Only a tiny portion of the parameter space, corresponding to a rather light SUSY 
spectrum and a conspiration of the IC^.J parameters, is already excluded. This result is 



consistent with the semi-analytical estimate provided in Eq. ( 111 ). The final sensitivity reach 
on BR(/j — > ej) of the MEG experiment [36] is represented by a black dashed line. We see 
that MEG is able to test indirectly a large part of the parameter space. As expected by the 



discussion in Sec. 4.5, BR(r — > ^7) cannot be as large as 0(10 ), thus resulting beyond the 
reach of the SuperB factory at KEK [62]. A positive signal for r — > fij would then disfavour 
our class of models. There is a corner of the parameter space which can be in the reach of a 
Super Flavour factory at the level of BR(r — > ^7) ~ 10 -9 [53]. Interestingly, these few points 
basically correspond to BR(/j — > &y) > 0(1O -13 ) and will be then tested by MEG, so that a 
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BR(x^|X y) 

Figure 8: BR(/i —y ej) versus BR(r — y fij). The thick black line represents the current 
experimental limit BR(/i —y 67) < 1.2 x 10 -11 [70], the dashed line the expected sensitivity of 
the MEG experiment. 



negative result would preclude the possibility of observing r — y fij at a Super Flavour factory 
as well. 

The fact that still a wide portion of the points in Fig. [8]would escape the searches at MEG 
is a consequence of our choice of scanning the parameter space up to very heavy slepton masses 
(~ 4 TeV). However, it is worth notice that even such a heavy spectrum can provide rates 
of (J, — y &y in the reach of MEG. Moreover, there are very good prospects to test most of 
the points of Fig. [8] in the future experiments looking for fj, — y e transition in nuclei. We 
remind that, within SUSY models, there is typically a striking correlation between /x — y e-y 
and the \x — y e conversion in nuclei, namely the fx — y e conversion rate is well approximated 
by CR(fi ->■ e in N) — Qiem x BR(/i — y 67)- This means that our prediction for BR(/i — y cy) 
can be easily translated in a prediction for CR(^x — y einN). The proposed fi— e conversion 
experiments at Fermilab [63] and at J- PARK [65J aim at the respective sensitivities of 10~ 16 
and 1CP 18 on CR(/x — y e in Ti). Therefore, such experiments would access the parameter 
space displayed in Fig. [8] almost completely, thus testing the model well beyond the MEG 
reach. 

As for LFV processes at colliders, such as X2 ~^ Xi/ xe ' we notice that they are very much 
suppressed in our scenario. This is due to the spectrum selected by the fermion masses 
and mixings fit, for which sleptons are heavy, and the second neutralino X2 mass is always 
smaller than the slepton ones, apart from some points with a light f±. Therefore, the LFV 
neutralino decay is a 3-body decay with off-shell intermediate sleptons, which give a GIM-like 
suppression of the process. 

In the left panel of Fig. |9j we plot BR(/i — y ej) versus the SUSY contribution to the 
anomalous magnetic moment of the muon a^ USY = (g — 2)^ USY /2. We see that no points 
provide a^ USY > 10~ 9 , a contribution which would lower the present tension between mea- 
surements and theoretical predictions below the 2a level. This is due to the heavy slepton 
spectrum favoured by the fit, which clearly determines a suppression of the SUSY contribu- 
tion to (g — 2)^. Within our model, a contribution at the level required by the present tension 
could be achieved in set-ups with larger tan /3, e.g. ~ 50 (since (g — 2)^ USY ~ tan 2 f3). 
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Figure 9: Left: BR(^ — > ey) versus the SUSY contribution to (g — 2) M /2. Right: BR(fi — > e^y) 
versus the electron EDM, d e . 



In the right panel of Fig. [9l we show an estimate for the electron EDM, d e . As discussed 



in Sec. 4.5 d e receives both flavour-dependent and flavour-independent contributions. Such 
an observable then results highly complementary to BR(/i — > &y) in testing the parameter 
space of the model, while it does not impose at present a further constraint to parameters 
of the model (the current limit being d e < 1.4 x 10 -27 e cm |66|). Given that the future 
experiments should reach d e ~ 10~ 30 e cm, or below [67], they could be able to test the 
model even beyond the reach of MEG, for vanishingly small values of BR(/i — > cy). On the 
other hand, a portion of the parameter space in the reach of MEG should escape the next 
generation eEDM experiments. However, it is interesting to notice that a large part of the 
parameter space should be in the reach both of MEG and the eEDM experiments. These 
two observables are therefore highly complementary and the interplay between both of them 
represents a powerful test of the class of models we are discussing. In fact, an evidence of 
H — > at MEG, at the level of, say, BR(/x — > ej) ~ 0(1O~ 12 ), would predict a lower bound 
on d e in the reach of the future experiments. On the other hand, observing the electron eEDM 
with a value of the order of 10 -29 e cm would imply that n — > e'y can be observed at MEG. 

Recent analyses of the consistency of the Unitarity Triangle (UT) have shown some ten- 
sions among the different observables used in the fit [681 169] . A way to solve such tension 
could be a positive new physics contribution to the Kaon CP-violating parameter ex at the 
level of about +20% [261 E3 EH]- As discussed in the previous sections, our class of models, 
being embedded in a SU(5) GUT, provides new sources of flavour violation both in the lep- 
tonic and the hadronic sector. In the left panel of Fig. 10, we show the predicted BR(/i — > ej) 
versus ex / ef^ • As we can see the SUSY contribution to ex can give sizeable deviations from 
the SM prediction. For instance, the cut of the left side of the plot reflects the fact that a 
large negative contribution to ex is excluded. A contribution of order +20% is still possible 
for a small portion of the parameter space. Moreover, most of such points, which correspond 
to values of BR(^ — > e-y) below the present limit, will be tested at MEG. 

B —7- tu, such as B s — > might represent an important constraint for SUSY models in 

the regime of medium or large tan j3. This is due to a potentially large charged Higgs mediated 
contribution, which grows with tan 4 f3 and always exhibits the opposite sign with respect to the 
SM contribution, whereas the present experimental measurement prefers values larger than 
the SM one (see, for instance, [26]). We plot BR(^ ey) and BR(B rv) /BR(5 tv) sm 



37 




Figure 10: fi — > versus the deviation from the SM predictions of ej< (left) and BR(B — > tv) 
(right). 



in the right panel of Fig. 10 We see that, even if most of the points give a prediction close to 
the SM one, it is still possible to have large deviations, which in future could be excluded by 
increasing the precision of the experimental determination of BR(i3 — > tv). Moreover, most of 
the depicted points lie in the 95% C.L. range 0.52 < BR(£ tv) /BR(£ tv) sm < 2.61 [26]. 
This is due to the Higgs spectrum selected by our model. In fact, m 2 H± (~ fj? + m 2 Hd + m^) 
cannot be too light since, as mentioned above, our large At gives sizeable \j? and tan/3 is 
not large enough to drive fn\j to large negative values. Only for few points of our sample, 
the parameters conspire to give a quite light Higgs sector. In fact, we find, as minimum 
values, rriA ~ 261 GeV, m H + ~ 272 GeV. These lower values correspond to points with 
quite large positive values of C^, which increase \Ab\ and \A T \ and therefore contribute to 
drive m# negative. These points give the largest deviation from the SM prediction for 
BR(£? — > tv) and some of them (cor responding to BR(£> tv)/BR(B -> tv) sm < 0.6) 
predict BR(1? S — > /x + /i~) at the level of 10 -8 , i.e. in the reach of the LHCb experiment in 
the near future. However, the bulk of the points just predicts no large deviation from the SM 
prediction, i.e. BR(B S /j, + /jl~) ~ O(10~ 9 ). 



6 Summary and Conclusions 

In this paper, we have investigated aspects of "SUSY flavour" models, towards predicting 
both flavour structures, in the context of SUGRA. We highlighted the importance of includ- 
ing carefully all the SUSY-specific effects such as one-loop SUSY threshold corrections and 
canonical normalization effects when fitting the model to the low energy data for the fermion 
masses and mixing angles. These effects entangle the flavour model with the SUSY parame- 
ters and leads to interesting predictions for the sparticle spectrum as well as, for instance, for 
the neutrino parameters. In addition, family symmetries introduced to explain the flavour 
structure of the Standard Model fermions can also make predictions testable in future flavour 
experiments. 

The extension of a pure model of flavour to a SUSY flavour model, and its phenomeno- 
logical analysis, may be performed in the following steps: 

• The starting point may be a family symmetry, as for instance SU(3), SO(3), or one of 
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their non-Abelian discrete subgroups like S4 or A4 (cf. Sec. [2]). 



The discussion of the SUSY formulation of the model starts by defining the superpo- 
tential W. 

The operators induced in the Kahler potential K by a given set of (flavour) messenger 
fields have to be included. 



• By canonically normalizing the fields (cf. Sec. 3.3) the canonically normalized Yukawa 
matrices are calculated. In this way, the structure of K enters the model predictions 
for the SM flavour parameters (i.e. the fermion masses and mixings). 

• From K and W, the (canonically normalized) soft SUSY-breaking parameters can be 
calculated. 

• To compare the flavour model predictions to the data, the RG evolution of all parameters 
has to be calculated. As discussed in Sec. 3.5., the inclusion of the one-loop SUSY 
threshold corrections is mandatory for all three generations at least for moderate and 
large tan f3 - and connects the flavour model analysis to the SUSY parameters. 

• Finally, the SUSY flavour structure, shaped by the same symmetries which govern the 
SM flavour structure, can give rise to various effects testable at flavour experiments. 

We have demonstrated these issues and applied the "SUSY extension" to a novel class of 
family symmetry models, assuming sequestering of W and K into visible and hidden sectors. 
The extension to SUSY, in addition to predictions for the SUSY flavour structure from the 
symmetries of the models, has various interesting effects: 



• Firstly, one-loop SUSY threshold corrections have turned out to have a crucial effect on 
the fit of the model to the data. Not only the ratios of quark and charged lepton masses, 
which are predicted by the model at the GUT scale but also the quark mixing angles, 
depend crucially on the threshold effects. Via them the SUSY spectrum enters into 
the fit of the model to the data and from this constraints on the SUSY parameters can 
be obtained. Using MC Monte Carlo techniques we demonstrated this for our example 
class of models. The main result is shown in Fig. [6j We also compared our full MC 
Monte Carlo analysis with a simple semi-analytical treatment, which may be applied as 
a first step to understand such kind of constraints on the SUSY spectrum from a SUSY 
flavour model analysis. 

• Secondly, also the canonical normalization of the kinetic terms, which generically re- 
quires a transformation (including rescaling) of the fields in a SUGRA context, modifies 
the Yukawa matrices and thus affects the fit. In our class of models the canonical nor- 
malization effects were dominated by a single parameter (related to a comparatively 
large ratio of the flavon vev ^3 over the messenger mass), a typical situation in flavour 
models with non-Abelian family symmetries (see, e.g. |14|). In our example class of 
models we saw, for example, that the prediction of tri-bimaximal mixing was modi- 
fied by this effect, which imposed constraints on the canonical normalization parameter 
C|/2, and predicts a correlation in the deviation from tri-bimaxmial mixing (cf. also 
|30| ) which could be tested in future precision neutrino oscillation experiments. 
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• Thirdly, the example class of flavour models we considered and generalized here turned 
out to be an interesting novel class of SUSY flavour models - since the flavour structure 
of the soft SUSY breaking terms turned out to realise a scheme we dubbed "Trilinear 
Dominance" . While the soft masses at the GUT scale have only very small off-diagonal 
elements, the flavour violating effects are mainly induced by the misalignment between 
the trilinear terms and the Yukawa matrices, governed by only a few parameters. 

We conclude that if SUSY will be discovered at the LHC, models of flavour would have to 
be extended to SUSY-flavour models. The testability of such models could be enhanced sig- 
nificantly, for example by the predictions for the sparticle spectrum and tan/3. Furthermore, 
flavour effects induced by the SUSY flavour structure could give new insights into flavour 
models, and characteristic SUSY flavour schemes (like "tri-linear dominance" discussed here) 
might be disentangled from others. Finally, even precision tests of neutrino mixing angles 
can be sensitive to SUSY/SUGRA typical features like canonical normalization corrections, 
as we have seen in our example class of models. 
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